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In  this  letter,  the  relationship  between  the  characteristic  function  for  two  arbitrary  noncommuting  observables  and  a 
generalized  Wigner  distribution  function  is  established.  This  distribution  function  is  shown  to  have  no  simple  interpretation 
in  the  sense  of  probability  theory  but,  in  lieu  of  its  special  properties,  can  be  used  directly  for  calculating  the  expectation 
values  of  observables. 


Whereas  classical  transport  physics  is  based  on  the 
concept  of  a  probability  distribution  function  which 
is  defined  over  the  phase  space  of  the  system,  in  the 
quantum  formulation  of  transport  physics,  the  con¬ 
cept  of  a  phase  space  distribution  function  is  not  pos¬ 
sible  inasmuch  as  the  noncommutation  of  the  posi¬ 
tion  and  momentum  operators  (the  Heisenberg  uncer¬ 
tainty  principle)  precludes  the  precise  specification  of 
a  point  in  phase  space.  However,  within  the  matrix 
formulation  of  quantum  mechanics,  it  is  possible  to 
construct  a  “probability”  density  matrix  which  is  of¬ 
ten  interpreted  as  the  analog  of  the  classical  distribu¬ 
tion  function. 

There  is  yet  another  approach  to  the  formulation 
of  quantum  transport,  based  on  the  construction  of 
the  Wigner  distribution  function  [1].  As  we  shall 
show,  this  distribution  function  has  no  simple  inter¬ 
pretation  in  the  sense  of  probability  theory  (2]  but, 
in  lieu  of  its  special  properties,  can  be  used  directly 
for  calculating  expectation  values  [3 — 5  J  of  observ- 

1  Partially  supported  by  the  Army  Research  Office. 
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ables  in  a  manner  quite  analogous  to  that  of  classical 
theory,  i.e.  by  integrating  the  product  of  the  observ¬ 
able  and  the  Wigner  distribution  function  over  all 
phase  space. 

This  letter,  in  part,  reviews  the  salient  features  of 
the  Wigner  distribution  function.  Although  the  Wigner 
function  is  generally  defined  in  terms  of  all  the  gener¬ 
alized  coordinates  and  momenta  of  the  system  in 
question  as 

/'w(xi...x„,p1...p„)  =  —  J  d/j...  dy„ 

+iyi,~.x„+\y„)  (0 

X  *(X,  -\ylt  ...x„  -\y„) 

X  exp[i(pjyj  + ...  +p„y„)/ft)  , 

we  will  discuss  the  properties  of  the  Wigner  function 
in  terms  of  a  single  coordinate  and  momentum.  In 
this  case,  we  let 
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iWC*.  p) -  ^  /  dy  *'(x+iy)*(x  -  {yypy'n 


(2) 


where  ♦(x)  refersto  the  state  of  the  system  in  the  co¬ 
ordinate  representation. 

The  distribution  function  of  eq.  (2)  has  interesting 
properties  in  that  the  integration  of  this  function  over 
all  momenta  leads  to  the  probability  density  in  real 
space;  conversely ,  the  integration  of  this  function  over 
all  coordinates  leads  to  the  probability  density  in  mo¬ 
mentum  space.  In  mathematical  terms, 


function  to  the  classical  distribution  function  is  appar¬ 
ent  by  examining  the  equation  of  time  evolution  for 
f*w(x,  p).  Upon  assuming  that  'I'fx)  in  eq.  (2)  satisfies 
the  Schrddinger  equation  for  a  system  with  hamilto- 
nian  H  -  p2/2m  +  K(x),  it  can  be  readily  shown  that 
/*w(x,  p)  satisfies  the  equation 

bPyjbt  +  (p/m)BPw/dx  *  0  •/>w  =  0,  (5) 

where 


8‘PW 


"  nm0 


(-iy 


(Jt/2)2"*1 
(2n+  1)! 


/  /,W(*.P)dP  =  ♦’(x^x) 


(3a) 


v32"*JK(x)  d2”"Pw(x,p) 
ax2"*1  dp2"* « 


(6) 


and 

/  Pv/(.x*P)  dx  -  0*(p)0(p) .  (3b) 

—  mm 

where 

«(p)  =  (2ir»)-1/2  J  e~iPJC/*^'(x)  dx  . 

—  mm 

It  follows  immediately  from  eq.  (3)  that,  for  an 
observable  W(x,p)  which  is  either  a  function  of  mo¬ 
mentum  operator  alone  or  of  position  operator  alone, 
or  any  additive  combination  therein,  the  expectation 
value  of  the  observable  is  given  by 


It  is  evident  that  in  the  limit  h  -*  0, 6  •  Py  in  eq.  (6) 
becomes 

«-Pw  =  -<d*73x)(3/>w/3p)  (7) 

so  that  eq.  (S)  reduces  to  the  classical  collisionless 
Boltzmann  equation. 

The  Wigner  distribution  function  defined  in  eq.  (2) 
is  derivable  (6)  from  the  Fourier  inversion  of  the  ex¬ 
pectation  value  (with  respect  to  state  ♦(x)  of  the  op¬ 
erator  e>(Tp+*x)  (here,  x  and  p  satisfy  the  commuta¬ 
tion  relation  [x,  p)  =  ih).  As  such, 

P)  =  ~// Cw(t,  •Je-KtP^dT  dO  ,  (8a) 

4rr^ 

where 


If  WPV/(x,  P )  dx  dp  , 


(4)  Cw(T,3)=yV(x)cM'P*®i>4'(x)dx,  (8b) 


which  is  analogous  to  the  classical  expression  for  the 
average  value.  Herein  lies  the  interesting  aspect  of  the 
Wigner  distribution  function;  the  result  of  eq.  (4) 
suggests  that  it  is  possible  to  transfer  many  of  the  re¬ 
sults  of  classical  transport  theory  into  quantum  trans¬ 
port  theory  by  simply  replacing  the  classical  distribu¬ 
tion  function  by  the  Wigner  distribution  function. 
However,  unlike  the  density  matrix,  the  Wigner  dis¬ 
tribution  function  itself  cannot  be  viewed  as  the  quan¬ 
tum  analog  of  the  classical  distribution  function  since 
it  is  generally  not  positive  definite  and  nonunique 
[P\f/{x,p)  of  eq.  (2)  is  not  the  only  bilinear  expres¬ 
sion  [1 ,3-5]  in  4'  that  satisfies  eq.  (3)] . 

Further  resemblance  of  the  Wigner  distribution 


and  the  interval  of  integration  is  (— °°,  °°)  unless  other¬ 
wise  specified.  In  order  to  show  that  the  right-hand 
side  of  eq.  (8a)  is  indeed  the  Wigner  distribution  func¬ 
tion  as  defined  in  eq.  (2),  note,  from  the  Baker- 
Hausdorff  theorem  [7],  that  ei(TP+e^*  can  be  written  a; 

ei(r p+ix)  -  eirp/2eiexeirp/2  f  (9) 

in  which  case  Cw(r,  6)  of  eq.  (8b)  becomes 


Cw(r,0)  =  f  [e-"P/2*(x))*e‘**[e-‘»P/2W(x)J  dx  , 
--  (>0) 


which  further  reduces  to 
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Cw(r,  e) =  /  **(x  -  i«*)c»*'»(*  +  \rh)  dx  .  (1 1) 

Then,  by  inserting  Cw(r,  6)  of  eq.  (11)  into  the  right- 
hand  side  of  eq.(8a),  integrating  over  the  variable  0 
by  using  the  relation 

f  d 0  =  2ir5(x'  -  x") , 

—  mm 

and  letting  r  =  -y/h,  the  desired  result  is  obtained. 

The  method  outlined  above  to  arrive  at  the  Wigner 
distribution  function  is  based  on  the  notion  of  a  char¬ 
acteristic  function.  The  characteristic  function  of  an 
observable,  A,  with  respect  to  state  140  (here,  the 
Dirac  notation  is  utilized  for  purposes  of  generality) 
is  defined  as 

C4tt)=<*lei{yil*>  ,  (12) 

where  { is  a  real  parameter.  Assuming  At  to  possess  an 
eigenvalue  spectrum  given  by  A  \A ')  =A'lA'),  CA($) 
can  be  evaluated  in  the  A  '-representation  as 

^(0=  fdA'f  dAK{^\A')<A'W^\A”)<A''m . 

(13) 

Since  C4'|cit<4  \A“)  5(4'  -4")  in  the  A '-rep¬ 

resentation,  C,|(€)  in  eq.  (13)  reduces  to 

CA(V  =  fdA'W  |^.|2,  (14) 

where  l^-l2  =  |<4'|40|2  =/*(4'),  the  probability 
distribution  function  for  measuring  A'  while  jn  state 
14').  Hence,  the  characteristic  function  for  A  is  the 
Fourier  transform  of  the  probability  distribution 
function  P(A’).  Subsequent  inversion  of  eq.  (14) 
above  leads  to 

/’(y,,)’^/C'«(t)e~iM  ‘ (,5) 

The  Wigner  distribution  function  was  derived  by 
taking  the  Fourier  transform  of  the  characteristic 
function  for  ei<TF+#*).  In  view  of  the  connection  be¬ 
tween  the  probability  distribution  function  and  the 
characteristic  function  for  a  given  observable,  this  ap¬ 
proach  seems  to  be  a  natural  way  of  obtaining  a  dis¬ 
tribution  function  for  momentum  and  position.  Un¬ 
fortunately,  the  noncommutative  nature  of  the  two 


observables  destroys  the  convenient  probability  inter¬ 
pretation  of  the  characteristic  function  implicit  in 
eq.(15). 

In  order  to  demonstrate  this  point,  assume  the  char¬ 
acteristic  function  for  two  noncommuting  observables, 
A  and  B,  to  be 

cA*(M2)=  .  (16) 

Observables  A  and  B  are  assumed  to  have  eigenvalue 
spectra 

A\A‘)-A'\A‘)  ,  B\B')^B'\B^) ,  (17) 

and  are  chosen  so  that  [A,  [A, B]  ]  =  [B,\A,B]]  - 0. 
This  assumption  is  imposed  so  that  the  identity 

e*<ti4+tafl)  =eit»'4*  eitl^e-Ut*M'^l/2  (18) 

may  be  used. 

Inserting  eq.  (18)  into  eq.  (16)  while  obtaining  the 
matrix  elements  of  ei*>>4  in  the  A  '-representation  and 
in  the  ^'-representation  results  in 

CAB(hh)  =  El  [AM  12 1  M'fdB’  (19) 

X  e^* a s ')  <4k  |y4  ')C4 '  |^’><^'  |  ^> . 

In  eq.  (19),  it  is  assumed  that  [A,B]  is  a  c-number  in¬ 
dependent  of  the  eigenvalues  A’  and  B'.  We  define 
F(A',  B'),  the  generalized  Wigner  distribution  func¬ 
tion*1,^  be 

F{A',B')=W\A'){A'\B'HB'\*),  (20) 

so  that 

F(A',B')=-L—  fdh  f  dbetihM*^)/2 
(2n)2  J  J 

XCABUlti2)'-iWWK  (21) 

It  is  evident  from  eqs.  (20,21)  that 


*’  The  form  of  the  generalized  distribution  function  derived 
is  sensitive  to  the  manner  in  which  exp|i({|4  +  tjA)J  is 
expanded.  l  or  example,  if  instead  of  eq.  (18),  we  used  the 
forms  exp(itjB) exp(it|A)exp(tifeM,  6)12),  exp0{,4/2) 
X  exp(it]B)exp(it)4/2),  or  exp(i{]£/2)exp(it|/4)  expfitj 
X  BIT),  all  of  which  are  equivalent,  we  would  indeed  ob¬ 
tain  a  different  form  of  the  generalized  distribution  func¬ 
tion,  yet  one  which  obeys  the  sum  rules  of  eq.  (22). 
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J>M\B')<M'  =  l(B'|*)|2  /d|2  (22a) 

X<Ua«U2)e-tt**' 

and 

(Mb) 


Thus,  eq.  (21)  establishes  the  relationship  between  the 
characteristic  function  for  two  arbitrary  noncommut¬ 
ing  observables  and  the  generalized  Wigner  distribution 
function.  The  generalized  distribution  function  has 
the  essential  properties  of  the  conventional  Wigner 
function  in  that  an  integration  of  the  generalized 
function  over  the  eigenvalue  spectrum  of  one  observ¬ 
able  leads  to  the  probability  density  in  the  canonically 
conjugate  observable  [eq.  (22)]. 

There  is  no  simple  probability  interpretation  of 
F(A',B?)  in  eqs.  (20, 21)  because  of  the  necessary 
overlap  between  the  states  of  the  noncommuting  ob¬ 
servables.  However,  if  A  and  B  are  made  to  commute 
so  that  \A')  and  IB')  are  a  common  set  of  eigenvec¬ 
tors,  then  F(A',B')  reduces  to  the  probability  distri¬ 
bution  function  for  A  and  B. 

Finally,  it  is  noted  that  the  conventional  Wigner 
distribution  function  for  observables  A  and  B  is  qb- 
tained  from 


(23) 

with  defined  in  eq.  (16),  whereas  the  alter¬ 

native  distribution  function, F(A',B’),  introduced  in 
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eqs.  (20, 21 )  differs  from  the  Wigner  function  due  to 
the  presence  of  the  phase  factor  e*i  fa  M  1/2  jn  the 
integrand  of  eq.  (21).  ?oiA=x  and  B  =p,Bw(x,p) 
in  eq.  (23)  reduces  to  the  Wigner  function  of  eq.  (2), 
whereas  F{x,p )  defined  from  eq.  (20)  becomes 

F{x,p)  =  ^  /  &  *•«* '{x-y)eWffl 


s(2  ff*)-I/2**(;r)eVJt/,,0(p),  (24> 

where  4>(p)  is  defined  in  eq.  (3b).  It  is  evident  that 
there  is  a  family  of  functions  which  are  bilinear  in  W 
yet  satisfy  the  sum  rules  of  eqs.  (3a,  b). 

There  are  some  interesting  questions  to  be  resolved 
concerning  the  uniqueness  and  positive  definiteness 
of  Wigner-type  quantum  distribution  functions.  Never¬ 
theless,  these  distribution  functions  serve  a  useful 
purpose  for  calculating  quantum  mechanical  observ¬ 
ables  in  transport  [S]  studies  and  numerous  solid- 
state  [8,9]  problems. 
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